Smectic Liquid Crystals in an Anisotropic Random Environment 
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We report a high-resolution x-ray scattering study of the smectic liquid crystal octylcyanobiphenyl 
(8CB) confined to ahgned colloidal aerosil gels. The aligned gels introduce orientational fields that 
promote long-range nematic order while imposing positional random fields that couple to the smectic 
density wave and disrupt the formation of an ordered smectic phase. At low densities of aerosil, 
the low-temperature scattering intensity is consistent with the presence of a topologically ordered 
XY Bragg glass phase that is predicted to form in response to such anisotropic quenched disorder. 
The observed features of the phase include an algebraic decay of the smectic correlations, which is 
truncated at large length scales due to the imperfect nematic order, and a power-law exponent that 
agrees closely with the universal value predicted for the XY Bragg glass. At higher aerosil densities, 
deviations from the XY Bragg glass form are apparent. At high temperature, the scattering intensity 
displays pre-transitional dynamic fluctuations associated with the destroyed nematic to smectic- 
A transition. The fluctuations obey quasi-critical behavior over an extended range of reduced 
temperature. The effective critical exponents for the correlation lengths and smectic susceptibility 
differ systematically from those of pure 8CB, indicating that coupling of the nematic order to the 
gel suppresses its role in the smectic critical behavior. 

PACS numbers: 61.30.Pq, 61.30.Eb, 64.70.Md, 61.10.Eq 



I. INTRODUCTION 

The effect of impurities and quenched disorder in con- 
densed matter represents an important problem both 
because of the ways in which disorder can change fun- 
damental properties of a system and because any real 
material will inevitably possess imperfections. Studies 
have shown dramatic consequences of disorder, such as 
the destruction of ordered phases and the introduction 
of new exotic ones, deviation from the universal behav- 
ior around transition points, and sometimes even the en- 
hancement of useful material properties. Liquid crystals, 
due to their soft elasticity, sensitivity to surface interac- 
tions, and experimental accessibility, provide an excellent 
opportunity for studying the effects of quenched disorder, 
which can be introduced experimentally through confine- 
ment of the liquid crystal to random porous media. A 
particular focus in this area has been the study of smec- 
tic liquid crystals in silica gels. In a detailed theoreti- 
cal treatment of smectics in such random environments, 
Radzihovsky, Toner, and coworkers identified two forms 
of random-field disorder that the confinement introduces 
to the smectic - random orientational fields that couple 
to the nematic director and random positional fields that 
couple to the smectic density wave [H, H, 0, |3) HI . One 
conclusion of their study was that the smectic phase is 
unstable to arbitrarily weak random-field disorder, con- 
sistent with expectations for a phase that breaks a con- 
tinuous symmetry in three dimensions. Experimentally, 
considerable work has been dedicated toward investigat- 
ing smectics in two types of silica gels, chemically rigid 
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and structures of the two gels are very similar, but the 
hydogen-bonded aerosil gels have the potential to access 
a weaker regime of disorder. Utilizing a wide range of 
probes, experiments on smectics confined to aerogels and 
aerosil gels have detailed the consequences of the confine- 
ment for the nematic and smectic phases, supporting the 
general picture that the gels impose random-field disor- 
der that destroys the nematic to smectic transition. 

Among the predictions for smectics confined to gels is 
the potential formation of topologically ordered, "Bragg 
glass" phases that are distinct from the high-temperature 
nematic phase by the absence of unbound dislocation 
loops. Considerable evidence from theory and simula- 
tion supports the existence of a topologically ordered 
state with algebraic decay of correlations in the three- 
dimensional (3D) random field XY model at low temper- 
ature [30, 3l|, [sl]. (However, a recent argument based on 
the functional renormalization group approach has ques- 
tioned the existence theoretically of such a phase in 3D 
XY systems [1^.) Since the nematic to smectic- A (N- 
SmA) transition breaks 3D XY symmetry, one might 
hence expect that smectics confined to gels would be 
good candidates for reafizing such an "XY Bragg glass" 
phase. Considering this possibility theoretically, Radzi- 
hovsky and Toner have concluded that the orientational 
random fields coupling to the nematic order and the soft 
elasticity of the smectic make smectics confined to gels 
distinct from standard random field XY systems. In- 
stead, they have introduced the possibility that a "smec- 
tic Bragg glass" phase, which has short range correlations 
and is qualitatively different from the XY Bragg glass, 
might form in smectics confined to gels However, 
the theoretical argument for the stability of the smectic 
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Bragg glass is inconclusive, and experimental evidence 
for the phase in smectics confined by silica gels has been 
contradictory @, [13, [13, [13] ■ Jacobsen, Saunders, et al. 
have further argued that smectics confined to uniaxially 
strained gels should fall into the universality class of ran- 
dom field XY systems and hence should form an "XY 
Bragg glass" for sufficiently weak disorder [3, [s^] ■ The 
key difference between isotropic gels and strained gels 
is the effect of the orientational fields. Specifically, if 
the strain alters the distribution of orientational fields, 
causing them to suppress nematic fluctuations at large 
length scales, then the smectic in strained gels should 
theoretically exhibit an XY Bragg glass, provided the 
positional random fields are sufficiently weak. Recent 
additional theoretical predictions regarding the smectic- 
A to smectic- C transition under confinement in strained 
gels has highlighted the potential for smectics in such an 
anisotropic environment to display exotic phase behav- 
ior m. 

In an effort to understand the separate effects of ori- 
entational random fields and positional random fields on 
the N-SmA transition and to test for the presence of an 
XY Bragg glass phase, we have conducted an x-ray scat- 
tering study of smectic octylcyanobiphenyl (8CB) con- 
fined to aligned aerosil gels. The structure of these gels 
dramatically alters the nature of the orientational fields, 
converting the random distribution found in isotropic gels 
into a sharply anisotropic distribution that macroscopi- 
cally orders the nematic director. This anisotropy, or 
"soft axis" , makes the gels a faithful experimental realiza- 
tion of uniaxially strained systems considered by Jacob- 
sen, Saunders, et al. In this paper, we report an analysis 
of the x-ray scattering results that indicates the presence 
of an XY Bragg glass at low temperature for low densi- 
ties of aligned aerosil. The signature feature of the XY 
Bragg glass is a power-law divergent smectic scattering 
peak with a universal exponent. We note that in a pre- 
liminary report on our x-ray studies we concluded 
that 8CB confined to aligned aerosil gels did not exhibit 
such characteristics of the XY Bragg glass phase. As we 
explain below, by adopting a physically better-motivated 
description of the behavior of the thermal critical fiuctu- 
ations in smectics confined to aerosil gels J28'], we find 
that the previous analysis is not satisfactory in account- 
ing for the smectic correlations at low temperature. In- 
stead, when a finite nematic domain size is incorporated 
into the analysis to refiect the imperfect nematic order in 
the aligned gels, comparisons between the x-ray scatter- 
ing results for low densities of aerosil and the predictions 
for an XY Bragg glass show excellent agreement. A dis- 
cussion of these findings is given in Sec. Ill below. 

Another important benefit of the macroscopic align- 
ment of the nematic director by the anisotropic gels is 
the pronounced enhancement in signal-to-background in 
x-ray measurements relative to measurements on smec- 
tics in isotropic gels. Taking advantage of this enhance- 
ment, we further provide in Sec. IV a detailed study of 
the smectic correlations at high temperature, which we 



analyze in terms of dynamic fluctuations that grow on 
approaching the critical region of the destroyed N-SmA 
transition. This analysis shows that the correlations, as 
characterized by effective critical exponents for the corre- 
lation lengths and susceptibility, differ from those of pure 
8CB even when the correlation lengths are considerably 
smaller than the gel pore size, providing evidence that 
the presence of quenched disorder alters the coupling be- 
tween nematic and smectic order and its influence over 
smectic critical behavior. When combined with similar 
results from smectics confined to isotropic gels for effec- 
tive critical exponents for the specific heat and smec- 
tic order parameter [l^, [H, [11], these findings provide 
a comprehensive picture of an evolution in the N-SmA 
critical fluctuations toward conventional 3D XY behav- 
ior with increasing quenched disorder. 



II. SAMPLE PREPARATION AND 
CHARACTERIZATION 

Pure 8CB undergoes an isotropic to nematic transi- 
tion at Tjyj = 313.98 K and a N-SmA transition at Tna 
= 306.97 K [l^. Confinement of the liquid crystal to 
aligned aerosil gels was accomplished through a proce- 
dure described previously [13]. Type 300 hydrophilic 
aerosil (DeGussa Corp.) was dried under vacuum at 393 
K for 24 hours. The dried powder along with appropri- 
ate quantities of 8CB (Frinton Laboratories, Inc.) was 
dissolved in high-purity acetone, and the mixtures were 
sonicated for at least two hours. The solutions were then 
heated to 318 K to evaporate the solvent slowly. After 
no visible trace of acetone remained, the samples were 
heated to 338 K under vacuum. The resulting composites 
consisted of 8CB confined to isotropic fractal aerosil gel s 
like those of previous studies 0, [ll [H, [H, [ll [H, [B] ■ 
To create anisotropic gels, we placed each sample in a 2 
Tesla magnetic field and cycled the temperature between 
308 K (nematic phase) and 318 K (isotropic phase) at 
least 100 times. Due to the magnetic anisotropy of ne- 
matic 8CB, the director tends to align parallel to the 
magnetic field, in competition with the random orienta- 
tional fields created by the gel. Due to the elasticity of 
the nematic and the compliance of the gel, this compe- 
tition results in restructuring of the gel to accommodate 
the magnetic anisotropy. As a result, the gels acquire a 
structure in which the orientational fields are no longer 
random but on average orient parallel to the aligning di- 
rection. Since the gel continues to have a random struc- 
ture positionally, the positional fields that couple to the 
smectic density wave remain random. Samples were pre- 
pared in this manner with aerosil densities ranging from 
Ps = 0.027 g sil/cni'^ 8CB, which is just above the gela- 
tion threshold, to ps = 0.10 g sil/cm^ 8CB, which is the 
upper limit above which efforts to align the gels were 
ineffectual. 

The smectic correlations that form in this anisotropic 
random environment were studied through high- 
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FIG. 1: Normalized scattering intensity of the smectic peak 
at fixed wave vector magnitude {q ~ 0.2 A~^) as a function of 
sample orientation for 8CB confined by an aligned aerosil gel 
with ps = 0.027 g/cm"^ (solid squares) and 0.10 g/cm"^ (open 
circles and solid triangles) at 298.15 K. The two scattering 
intensities for ps = 0.10 g/cm^ correspond to measurements 
before (open circles) and after (solid triangles) annealing the 
sample in the isotropic phase of 8CB. 



resolution x-ray scattering. The experiments were con- 
ducted at the X22A beam line of the National Syn- 
chrotron Light Source. The beam line is equipped with 
a Ge(lll) monochromator to select a beam of 10 keV x- 
rays. The measurements were performed in transmission 
geometry. The beam size was approximately 1 mm x 
1 mm, and the sample thickness was approximately 1.5 
mm to match the attenuation length of the 10 keV x-rays. 
A triple-bounce Si(lll) channel-cut analyzer crystal was 
positioned between the sample and a scintillation point 
detector to achieve high wave-vector resolution. 

To determine the degree of anisotropy in the orienta- 
tional fields coupling to the liquid crystal, we character- 
ized the mosaic spread of the smectic layer normal by 
"rocking curve" measurements through the smectic scat- 
tering peak. Figure [T] shows mosaic scans for ps = 0.027 
and 0.10 g/cm^ at T=:298.1 K, several degrees below 
Tna- The smectic scattering displays strong azimuthal 
anisotropy, indicating a narrow distribution of smectic 
layer-normal orientations. The alignment is slightly bet- 
ter (the peak is narrower) for lower density gels, consis- 
tent with the gels' increased ability to resist restructur- 
ing with increasing density due to greater in mechani- 
cal strength [l^, [2^ . We note these measurements were 
performed several weeks after removal of the magnetic 
field used to align the gels; thus, the observed macro- 
scopic alignment is a consequence of the gel and not 
the external field. As also shown in Fig. [H the qual- 
ity of the nematic alignment displays negligible change 
in response to temperature excursions into the isotropic 
phase, demonstrating that the anisotropic gel structure is 
stable against thermal fiuctuations and sufficiently robust 
to re-align the nematic director on subsequent cooling. 
Indeed, samples are observed to remain aligned after re- 
peated heating into the isotropic phase and after storage 
for several months. Thus, the random orientational fields 



of the isotropic gel have been converted into fields with 
long-range order that align the nematic director macro- 
scopically. 



III. SMECTIC CORRELATIONS: EVIDENCE 
FOR AN XY BRAGG GLASS 

Figure 2(a) displays the low-temperature x-ray scatter- 
ing intensity for ps = 0.027 g/cm^ as a function of wave 
vector parallel to the alignment direction [i.e., the direc- 
tion of the magnetic field used to prepare the gel), g^, 
and perpendicular to it, q^. In the alignment direction 
the intensity displays a pronounced peak corresponding 
to the smectic scattering, while in the perpendicular di- 
rection the intensity shows no evidence of smectic scat- 
tering, consistent with the strong anisotropy implied by 
Fig. 1. We ascribe the scattering intensity along en- 
tirely to the gel structure and note that an essentially 
identical background intensity contributes to the scat- 
tering along Qz . In order to isolate the smectic scattering 
from this background, we subtract the intensity along 
Qx from the intensity along qz. This procedure assumes 
that the anisotropy in the gel structure does not affect 
the scattering from the gel in the wave-vector range of 
the smectic peak. As described previously, small angle 
scattering measurements support this assumption [ij]. 
Nevertheless, the observed smectic signal greatly exceeds 
this background contribution, by a factor of more than 
a factor of 10"* in Fig. 2(a), and this strong signal-to- 
background, which results from the azimuthal focusing 
by the aligned orientational fields, enables the detailed 
examination of the smectic correlations described below. 
Figure 2(b) displays the results for the scattering inten- 
sity along Qz from Fig. 2(a) with the background sub- 
tracted. Also in the figure is the resolution function de- 
termined from the profile of the direct x-ray beam. The 
inset to Fig. 2(b) compares the scattering to the resolu- 
tion in the region of the peak. The peak is clearly broader 
than the resolution even for this lowest aerosil density at 
low temperature. 

As a first effort to characterize the smectic order, we 
consider a correlation function based on the random- 
field model that has been used previously to describe the 
short-range smectic correlations observed under confine- 
ment in isotropic aerosil gels [l3, [13, HE 113 ■ Specifically, 
we model the smectic peak with the two-component line 
shape 
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where gy and q± are the wave vectors parallel and per- 
pendicular to the smectic layer normal, respectively, and 
go is the ordering wave vector. The first term in /(q) 
is an anisotropic Lorentzian with fourth-order correction 
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FIG. 2: (a) X-ray scattering intensity for 8CB confined by 
an aligned aerosil gel with ps = 0.027 g/cm^ at 296.15 K 
as a function of wave vector parallel to the gel's alignment 
direction, (circles) and perpendicular to this direction, qx 
(triangles), (b) Scattering intensity along q^ with the inten- 
sity along qx subtracted to isolate the smectic scattering. The 
solid line in (b) is the instrumental resolution obtained from 
the incident beam profile. The inset to (b) shows the peak 
region on an expanded scale. 



that describes critical dynamic fluctuations on approach- 
ing the N-SmA transition in pure liquid crystals, ai is 
the smectic susceptibility, and ^ly and ^ij^ are the corre- 
lation lengths of the fluctuations in the directions parallel 
and perpendicular to the nematic director, respectively. 
The second term, which has a shape that is proportional 
to the square of the thermal fluctuation term, describes 
static short-range fluctuations due to the random fields, 
which are characterized by the correlation lengths ^2|| a-nd 
^2±. Previous studies of smectics in isotropic aerosil gels 
have found that the correlations display two temperature 
regimes. At low temperature, a2 > 0, and the amplitude 
of the static fluctuation term rises from zero with de- 
creasing temperature. At higher temperature, 02 — 0, 
and the scattering proflle is similar to that of the ne- 
matic phase of pure 8CB, where pre-transitional smectic 
critical fluctuations dominate O [13, HI] . 

Figure [3] shows the results of fits to Eq. ([T]) with 02 = 
at two temperatures in this high-temperature regime for 
8CB confined by an aligned gel with ps = 0.10 g/cm'^. To 
perform fits to Eq. ([1]), we integrate numerically over the 
distribution of layer normal orientations, as determined 
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FIG. 3: Background-subtracted scattering intensity of the 
smectic peak for 8CB confined to an aligned aerosil gel with 
Ps = 0.10 g/cm^ at (a) 309 K and (b) 306 K. The solid lines 
are the results of fits to Eq. ^ with 02 = 0. 



by rocking curves like those shown in Fig. [1] and convolve 
with the instrumental resolution. In this analysis and 
^ij^ are treated as independent parameters, while the am- 
plitude of the fourth-order correction term, c, is treated 
as a function of 1 with c(fix) set by the behavior in 



pure 8CB 



(Allowing c to vary as an addi- 



tional free parameter leads to significant scatter in the 
fit parameters.) As Fig. [3] illustrates, the dynamic term 
in Eq. ([1]) provides an excellent description of pre- tran- 
sitional smectic fluctuations in 8CB conflned to aligned 
aerosil gels at high temperatures. This agreement al- 
lows us to study in detail the pseudo-critical behavior of 
8CB with anisotropic quenched disorder, as described in 
Sec. IV below. 

Turning to low temperature, the dynamic fluctua- 
tion term in Eq. (1) alone cannot describe the mea- 
sured lineshape. As mentioned above, previous anal- 
ysis of the smectic correlations under confinement in 
isotropic aerosil gels [13, El, [H, H^l have modeled the 
low-temperature scattering by including the contribu- 
tion in Eq. (1) from static fiuctuations (02 > 0). In 
most cases, this analysis has assumed that ^2|| = Ci|| 
and ^2± = Ciij that is, the correlation lengths in the 
static and dynamic terms were assumed equal. Typi- 
cally, this approach provided good agreement with the 
measured line shapes. Indeed, in a preliminary report 
on 8CB confined to aligned aerosil gels, we found that 
Eq. ^ with f2|| = and ^2J. = £,i± described the data 
accurately [14| . However, this agreement does not neces- 
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sarily indicate that the analysis is capturing the correct 
physical picture. As pointed out recently by Larochelle 
et ai, the results of such analysis lead to trends in the 
smectic susceptibility and correlation lengths that are in- 
consistent with_calorimetric studies of smectics confined 
to aerosil gels [28|. Calorimetry measurements on 8CB 
in isotropic aerosil gels at low density {ps < 0.10 g/cm'^) 
show a sharp peak in the heat capacity, with power-law 
temperature dependence resembling pre-transitional crit- 
ical behavior extending over an extended range of re- 
duced temperature both above and below an effective 
transition temperature T* (l6j . This behavior would sug- 
gest that ^i_L, and cti should similarly decay as a 
power law as a function of reduced temperature away 
from T*. As discussed in Sec. IV, for T > T* the smectic 
fluctuations in 8CB confined to aligned aerosil gels indeed 
show a quasi-critical increase toward T*. However, anal- 
ysis of the x-ray results at low temperature using Eq. ([T]) 
with ^2|| — Ci|| and ^2± = Ci± leads to correlation length 
and susceptibility values that remain large away from T* . 
For example, ^|| « 3000 A at T*-T = 9 K for 8CB con- 
fined to an ahgned gel with ps — 0.10 g/cm^ [l3|, incon- 
sistent with pseudo-critical behavior. Such large dynamic 
fluctuations persisting to low temperature seem very un- 
likely. However, if we restrict the dynamic fluctuations to 
follow a more plausible temperature dependence, we find 
that Eq. (1) cannot describe the x-ray lineshapes for 8CB 
in aligned gels at low temperature. Rather, as we explain 
below, the lineshape predicted for the XY Bragg glass 
with an appropriate cut-off for a finite nematic domain 
size can very accurately describe the measured lineshapes 
for low densities of aligned aerosil. 

To examine the situation more closely, we display in 
Figs. (Ha) andlDJb) log- log plots of the scattering intensity 
versus {qz — qo) at various temperatures for ps = 0.042 
g/cm'^. As described in Sec. IV, the effective transition 
temperature for = 0.042 g/cm^ is T* = 305.0 K. As 
seen in Fig. |3]Ja), the intensities at temperatures above 
T* overlap at large {qz — qo)- This collapse indicates 
that the intensities away from the peak share the same 
wave-vector dependence, specifically that of the critical 
fluctuation term in Eq. ([T]), and further that the ratio 
5^11 /(Ti is a constant independent of temperature (imply- 
ing the critical exponent {2 — rj) w 2). Similarly, for the 
two temperatures just below T* shown in Fig. Ufa), 303.9 
K and 304.6 K, the intensities at large {q^ — qo) collapse 
onto the data for T > T*. This collapse indicates that 
critical fluctuations continue to dominate the intensity 
at large {qz — qo) at temperatures near but below T*. 
This finding is not surprising since one would expect to 
have large critical fluctuations in the vicinity of a true 
critical point. However, for lower temperature, the situ- 
ation becomes more complicated. Plotted in Fig. [H^b) is 
the intensity at 300.1 K, approximately 4.9 K below T*, 
along with the intensity at 309.1 K, approximately 4.1 
K above T*. Not only does the low-temperature curve 
fail to collapse onto the high-temperature curve at large 
{qz — qo), but it exhibits considerable excess intensity at 
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FIG. 4: Background-subtracted scattering intensity at various 
temperatures as a function of the difference in wave vector 
from the ordering wave vector, qo, for ps — 0.042 g/cm^. (a) 
intensities at several temperatures ranging from just below T* 
= 305.0 K to several degrees above T*. (b) intensities at two 
temperatures, 300.1 K and 309.1 K, that are approximately 
the same difference in temperature from T* . The dashed line 
in (b) is the result of a fitting the intensity at 300.1 K in the 
peak region {{q^ — qo) < 10~^ ^~^) to the second term in 
Eq. IT]) (di — 0) representing static short-range order. 



smaller {qz — qo) that cannot be accounted for by either 
term in Eq. ([1]). Specifically, the dashed line in Fig. UJb) 
shows the result of fitting the low-temperature intensity 
in the peak region to the static term of Eq. ([T|) alone 
((Ti = 0). While the fit agrees adequately with the mea- 
sured intensity in the immediate vicinity of the peak, it 
deviates significantly at larger {qz — qo)- A comparison 
with the high-temperature intensity demonstrates that 
the dynamic term in Eq. ([T]) cannot account for these 
deviations. Since the two intensities plotted in Fig. IH^b) 
are for temperatures that are approximately the same dif- 
ference in temperature from T*, they correspond to ap- 
proximately the same reduced temperature. Hence, the 
pre-transitional critical fluctuations at the two temper- 
atures should be similar (assuming the amplitude ratios 
for the correlation lengths and susceptibility are of or- 
der one). For scattering from such critical fluctuations, 
a power law in the intensity should be observed only for 
{qz — qo) > An inspection of the line shape for 

309.1 K indicates that for this reduced temperature this 
asymptotic range is reached only for {qz — go) > 4 x 10~^ 
A~^. The extension to lower {qz — qo) of the power-law 
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behavior at 300.1 K thus indicates that the scattering in- 
tensity is inconsistent with the small dynamic correlation 
lengths expected at a temperature so far below T*. In 
essence, the constraint used in the previous analysis 
that ^2|| = Ci|| s-iid ^2J_ = £,i± led to artificially large val- 
ues of and at low temperature that fortuitously 
approximated the measured line shape but that on fur- 
ther inspection cannot be justified as credible behavior 
for critical fluctuations. 

The extended power-law tails in the low-temperature 
scattering peak instead suggest the possible presence of 
the XY Bragg glass phase. As mentioned above, the XY 
Bragg glass is ideally characterized by a power- law diver- 
gent smectic scattering peak with a universal power-law 
exponent. As noted by the inset to Fig. 2(a), the mea- 
sured scattering peaks are appreciably broader than the 
instrumental resolution, which precludes such power-law 
divergence extending to the lowest {qz — qo). However, 
we argue that this discrepancy could be a consequence of 
a large length scale cut-off introduced by the liquid crys- 
tal ordering in the aligned gels. Specifically, while the 
mosaic scans in Fig. 1 indicate long-range nematic or- 
der, they also show an appreciable spread in the smectic 
layer-normal orientations about the aligning direction. 
This spread indicates spatial variations in the smectic 
layer normal that are likely set by the length scale char- 
acterizing the nematic order. In other words, the liquid 
crystal in the aligned gel has finite-size nematic domains 
like in the isotropic gels [s^l, the difference being that 
the domains are not randomly oriented but rather have 
a distribution of orientations clustered around the align- 
ing direction. Nevertheless, smectic order cannot persist 
across two misaligned domains, leading directly to a cut- 
off in the smectic correlations at low (g^ — qo). 

Taking such an effect into account, we model the scat- 
tering peak at low temperature with a form 



/(q) = P(q) + Hiq) 



(2) 



where P(q) is the XY Bragg glass correlation func- 
tion truncated at small (q^ — Qo), and H{q) accounts 
for the finite-size effect within a Gaussian approxima- 
tion [H, HOl . The XY Bragg glass form is given by [1, [s^l 
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if ((911 -qo)^+aql) > 1/L, 



(3) 



if ((911 ~qo)^+aql) < l/L. 



where 5 is a universal, temperature-independent expo- 
nent predicted to be 5 = 2.45 d, [s^l, C sets the over- 
all amplitude, a is a non-universal constant that is ex- 
pected to be of order unity [s^l , and L is the size scale of 
the nematic domains. The finite-size correction has the 
form [4(1 



i?(q) = Ae--^^(('?i|-9o)^+"9i)/4'r 



(4) 




where the amplitude A is related to C such that at wave 
vectors {{q\\ — qo)^ + aq\) < 1/L the Gaussian form of 
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FIG. 5: Comparison between the background-subtracted scat- 
tering intensity as a function of wave vector for ps = 0.027 
g/cm^ at 296.1 K and the form predicted for an XY Bragg 
glass with a truncation in the correlations at large length 
scales. The circles are the measured intensity, and the solid 
line is the result of a fit to Eq. (2). The inset shows a com- 
parison between the scattering intensity and the fit result in 
the peak region on an expanded scale. The XY Bragg glass 
form agrees closely with the measured intensity over the full 
range of wave vectors. 



H{q) dictates the line shape, while for larger wave vec- 
tors it has negligible effect on the line shape [s^. To 
perform fits to Eq. (2), we again integrate numerically 
over the distribution of layer normal orientations, as de- 
termined by rocking curves like those shown in Fig. [Jl 
and convolve with the instrumental resolution. Figure 5 
shows the result of a fit to Eq. (2) for ps — 0.027 g/cm'^ 
at T = 296.1 K, the lowest measurement temperature. 
The figure displays the scattering intensity over the full 
range of wave vectors, while the inset shows the peak re- 
gion on an expanded scale. The agreement between the 
measured line shape and the XY Bragg glass form is es- 
sentially perfect. Good agreement is similarly found for 
Ps = 0.42 g/cm^^, as illustrated by Fig. 6 which shows the 
result of a fit using Eq. (2) to the scattering intensity at 
T = 298.1 K for this aerosil density. The fits for both 
densities give L sa 1 pm, in reasonable agreement with a 
value expected based on the nematic correlation lengths 
measured in isotropic aerosil gels [s^. The anisotropy 
factor is found to be a w 0.15 in the two cases. For larger 
densities of aligned aerosil, the XY Bragg glass form is 
less successful in describing the measured line shape, as 
illustrated by Fig. 7 which shows the result of a fit using 
Eq. (2) to the low- temperature scattering intensity for 
Ps = 0.10 g/cm^. We interpret this discrepancy as an 
indication that the larger aerosil densities impose posi- 
tional disorder that is too strong for the XY Bragg glass 
to be stable. 

For the lower aerosil densities, the close agreement with 
the XY Bragg glass form extends over a range of temper- 
atures below T*. However, at temperatures sufficiently 
close to T* , contributions to the scattering from smectic 
critical fluctuations become appreciable. These contribu- 
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FIG. 6: Comparison between the background-subtracted scat- 
tering intensity as a function of wave vector for ps — 0.042 
g/cm^ at 298.1 K and the form predicted for an XY Bragg 
glass with a truncation in the correlations at large length 
scales. The circles are the measured intensity, and the solid 
line is the result of a fit to Eq. (2). The inset shows a com- 
parison between the scattering intensity and the fit result in 
the peak region on an expanded scale. The XY Bragg glass 
form agrees closely with the measured intensity over the full 
range of wave vectors. 
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FIG. 7: Comparison between the background-subtracted scat- 
tering intensity as a function of wave vector for pa = 0.10 
g/cm^ at 298.1 K and the form predicted for an XY Bragg 
glass with a truncation in the correlations at large length 
scales. The circles are the measured intensity, and the solid 
line is the result of a fit to Eq. (2). Deviations between the 
XY Bragg glass form and the measured intensity are apparent 
at large wave vectors. 



tions are apparent in the two temperatures slightly below 
T* included in Fig. 4(a) where both smectic critical fluc- 
tuations at large (q^ — go) E^nd additional contributions 
at smaller wave vector are apparent. Efforts to account 
for the contributions from dynamic fluctuations in this 
temperature region by adding a term to Eq. (2) are suc- 
cessful in the sense that good agreement with the mea- 
sured scattering intensity can be achieved. However, this 
analysis leads to large uncertainties in the fit parameters; 
therefore, we will not discuss this temperature region and 
will restrict our comparison to the XY Bragg glass form 



FIG. 8; Power-law exponent S for the XY Bragg glass as a 
function of temperature for 8CB confined to ahgned aerosil 
gels with densities ps = 0.027 g/cm^ (sohd) and 0.042 g/cm^ 
(open). The dashed line is the universal value for the expo- 
nent, S = 2.45, predicted theoretically [5ll34|. 



to lower temperatures, where the dynamic fluctuations 
can be safely neglected. Figure 8 displays the values of 
the power-law exponent i5 extracted from fits over this 
low-temperature region for ps = 0.027 and 0.042 g/cm^. 
The measured exponents for both densities are remark- 
ably close to the universal value predicted for the XY 
Bragg glass, S — 2.45 [1, [s^l, shown by the dashed line in 
the figure. We note that like the XY Bragg glass phase, 
the smectic phase of pure liquid crystals possesses quasi- 
long range order due to the Landau-Peierls instability 
and hence also displays a power-law singular scattering 
peak. However, the power-law exponent along gn for the 
smectic phase is always less than 2 [i^ |4l[ , a value that 
is incompatible with the measured line shapes for 8CB 
confined to aligned aerosil gels. In contrast, the close 
quantitative agreement between the exponent predicted 
for the XY Bragg glass and the observed values lends 
strong support to the conclusion that this topologically 
ordered phase forms in 8CB confined to aligned aerosil 
gels with low density. 

Despite this excellent, quantitative agreement between 
the measured scattering profiles for low densities of 
aerosil and the XY Bragg glass form, we note that some 
caution is in order. As described previously, one can ana- 
lyze the scattering intensity using Eq. (1) with ^2|| = ^i|| 
and £,2± — £,i±- (Although, the agreement is not nearly 
as accurate as that of the XY Bragg glass form nor does 
it cover as large a wave- vector range; see Fig. 5(c) in 
Ref. [l3|)- Further, the values that are obtained for 
^11 and from fitting to Eq. (1) are similar to those 
found for 8CB confined to isotropic aerosil gels follow- 
ing the same analysis As we argue above, this 
analysis provides a flawed physical picture; however, this 
coincidence in the correlation lengths suggests that the 
smectic correlations that form in the aligned gels and 
the isotropic gels are similar. Since the XY Bragg glass 
phase is not expected to be relevant to smectics confined 
to isotropic gels, one might interpret this similarity as 
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evidence against the XY Bragg glass phase forming in 
the aligned gels. Additional experimental work on other 
smectic liquid crystals confined to aligned gels that test 
for the presence of the XY Bragg glass phase would help 
in clarifying the stability of this topologically ordered 
phase within the anisotropic random environment cre- 
ated by the gel. 

IV. HIGH-TEMPERATURE CORRELATIONS: 
PSEUDO-CRITICAL BEHAVIOR 

As described above, the high-temperature scattering 
intensity of 8CB confined to aligned aerosil gels is char- 
acterized by dynamic critical fluctuations. We interpret 
these fluctuations as the remnants of the N-SmA critical 
point that is destroyed by the quenched disorder (and 
not directly with the transition to the XY Bragg glass 
phase which presumably has a significantly subtler sig- 
nature in the x-ray scattering than the observed robust 
peaks illustrated in Fig. 3). Although the N-SmA criti- 
cal point is obtained strictly only in the presence of zero 
quenched disorder, upon cooling in weak disorder the sys- 
tem in some sense comes in close proximity to the critical 
point and should therefore display critical smectic fluc- 
tuations. Further, if the disorder is sufficiently weak, the 
susceptibility and correlation lengths charactering these 
fluctuations, while non-singular, should display pseudo- 
critical behavior over an extended range of reduced tem- 
perature. This premise leads us to define the effective 
transition temperature T* as the temperature at which 
these quantities would diverge if their growth were not 
truncated: 

qo^m - ir"" (5) 

9o6± - \tr^ (6) 

where t = (T— T*)/T* is the reduced temperature, and 
v^, i'±, and 7 are effective critical exponents for ^ij^, 
and (Ti, respectively. 

If this assumption of truncated power-law divergences 
is valid, then the three critical parameters should yield 
self-consistent values for T* . To test this assumption and 
to identify the correct value of T* , we fit the values of || , 
and cTi to Eqs. I©-©. In these fits, only data at high 
temperatures where the dynamic term in Eq. (1) alone 
describes accurately the measured line shapes were in- 
cluded. As mentioned above, at lower temperature (near 
but below T*), additional contributions to the scatter- 
ing intensity, presumably from the static correlations of 
the XY Bragg glass phase, appear. Since efforts to sep- 
arate these different contributions lead to considerable 
uncertainties in the fit parameters, we focus on the high 
temperature side of the critical region. In the fitting, the 
critical exponents were held fixed and T* was treated as 
a free parameter. These fits were repeated over a range 
of critical exponents, and the quality of the fits, as deter- 
mined from the values, was compared. For all aerosil 



TABLE I: Summary of effective critical exponents for 8CB 
confined to aligned aerosil gels in comparison with pure 8CB 
and the 3D XY model. Also included are the gel pore chords 
lo Jj^ and effective transition temperatures T*. 



System 


lo (A) 


T* (K) 


7 


^11 






Pure 8CB 




306.97 


1.26 


0.67 


0.51 


1.314 


0.027 g/cm 


^ 3000 


305.96 


1.55 
±0.05 


0.73 

±0.02 


0.70 
±0.02 


1.05 


0.042 g/cm^ 


~ 1500 


305.01 


1.54 
±0.05 


0.74 
±0.03 


0.65 
±0.03 


1.13 


0.071 g/cm^ 


~ 900 


305.62 


1.56 
±0.05 


0.75 
±0.02 


0.66 
±0.02 


1.15 


0.10 g/cm^ 


~ 600 


305.86 


1.53 
±0.05 


0.73 
±0.02 


0.65 
±0.03 


1.13 


3D XY 






1.32 


0.67 


0.67 


1 



densities, the minima in for the three power-law rela- 
tions, Eqs. ([5])- ([7]), occur at values of T* that are within 
0.01 K of each other. We equate T* with the average of 
the three values determined in this way for each aerosil 
density. This value is listed in Table HI Comparing T* 
for different aerosil densities, we note the value displays a 
non-monotonic dependence on ps with a minimum near 
Ps = 0.05 g/cm"^, consistent with the effective critical 
temperature determined with calorimetry for 8CB con- 
fined to isotropic aerosil gels [l^ . 

Figures [9] and [To] display log - log plots of ^ly, and 
CTi as a function of reduced temperature for T > T* for 
various gel densities using the optimal values of T*. The 
solid lines in the figures are the results of fits to Eqs. ([5])- 
([7|) , which demonstrate that the temperature dependence 
is well described by power laws, further reinforcing the 
idea of pseudo-critical behavior. The effective critical 
exponents extracted from the fits are listed in Table |I] 
along with the values for pure 8CB [HllEIll]. The val- 
ues of i'±, and 7 for 8CB confined to aligned aerosil 
gels are distinctly different from those of pure 8CB, indi- 
cating that the quenched disorder introduced by the gel 
not only affects the low-temperature smectic order but 
also modifies the pre-transitional critical fluctuations. At 
sufficiently high temperature the correlation lengths are 
small compared to the mean pore chord Iq of the gels [l^ , 
which is also listed in Table U Hence, the disorder caused 
by the gels should have little direct effect on the smectic 
fluctuations well above T*. In the simplest picture of the 
N-SmA transition, the formation of the smectic density 
wave breaks 3D XY symmetry, and thus the critical be- 
havior of pure liquid crystals could be expected to match 
that of the 3D XY model. However, the observed critical 
behavior of pure 8CB and other smectic liquid crystals 
typically deviates from that of the 3D XY model. These 
deviations are illustrated by the difference in critical ex- 
ponent values for pure 8CB and the 3D XY model listed 
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FIG. 9; Dynamic correlation lengths ^1 II (circles) and ^ix (tri- 
angles) as a function of reduced temperature for 8CB confined 
to aligned aerosil gels with densities (a) ps — 0.027 g/cm^, (b) 
Ps = 0.042 g/cm^ (c) Ps = 0.071 g/cm^ and (d) ps = 0.10 
g/cm^. The solid line are the results of power- law fits. 
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FIG. 10: Smectic susceptibilty as a function of reduced tem- 
perature for 8CB confined to aligned aerosil gels of density for 
(a) Ps = 0.027 g/cm^ (b) ps = 0.042 g/cm^ (c) ps = 0.071 
g/cm^, and (d) ps — 0.10 g/cm^. The solid line are the results 
of power-law fits. 



in Table m The precise nature of these deviations, whose 
source likely involves couplings between the smectic or- 
der parameter and the nematic order, has been a long- 
standing problem. As described above, when a smectic 
liquid crystal is confined to an aligned aerosil gel, it ex- 
periences not only positional disorder that couples to the 
smectic density wave but also anisotropic orientational 
fields that couple to the nematic order, thereby modify- 



ing the nematic behavior. We therefore conclude that 
differences in the critical fluctuations between pure 8CB 
and 8CB confined to aerosil gels originate from the cou- 
pling of nematic order to the gel and its indirect influence 
on the smectic critical behavior. 

As an illustration, compare the ratio i^\\/i'± for vary- 
ing aerosil densities with the ratio for pure 8CB listed in 
Table [D For the 3D XY model, the ratio is one due to 
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FIG. 11: Critical exponents (a) (b) and (c) 7 for dif- 
ferent smectic liquid crystals as a function of McMillan ratio, 
Rm = • The values for various pure liquid crystals, shown 
by circles, are taken from [4§]. The values for pure 8CB are 
shown by the solid circles. Also included are the values of the 
effective critical exponents for 8CB confined to aligned aerosil 
gels with densities ps — 0.027 g/cm^, (triangle), ps — 0.042 
g/cm'' (diamond), ps — 0.071 g/cm'^ (upside-down triangle), 
and ps = 0.10 g/cm^ (square), which are assigned effective 
McMillan ratios using Eq. The solid lines mark the ex- 
ponent values for the 3D XY model, and the dashed lines are 
guides to the eye. 



the isotropic scaling of correlations. However, for pure 
8CB the scalings are highly anisotropic (j^||/i^_L = 1-314), 
consistent with behavior seen in various other smectic liq- 
uid crystals Patton and Andereck have advanced a 
theory that this anisotropic scaling results from coupling 
between nematic director fluctuations and the smectic 
order parameter [H, 13] ■ For 8CB confined to ahgned 
aerosil gels, the ratio i'\\/t^± shows much weaker scaling 
anisotropy than in pure 8CB. One interpretation for this 
decreased scaling anisotropy is that pinning of the ne- 
matic director by the anisotropic orientational fields sup- 
presses long-wavelength nematic fluctuations, thus weak- 



ening their influence over the smectic critical behavior. 

While the decrease in i'\\/i'± for 8CB confined to 
ahgned aerosil gels compared with i^\\/i^± for pure 8CB 
indicates qualitatively that the quenched disorder weak- 
ens nematic-smectic coupling, this observation would 
be made stronger with a quantitative measure of this 
weakening. In pure liquid crystals, the strength of the 
nematic-smectic coupling can be crudely parameterized 
by the McMillan ratio 



Rm 



Tna 
Tni 



(8) 



specifying the temperature range of the nematic phase. 
Large Rm indicates a short nematic range, which typi- 
cally implies the nematic order is far from saturated at 
the N-SmA transition. Hence, for large Rm the nematic 
susceptibility can be expected to be large, and concomi- 
tantly the order parameter coupling can be expected to 
be strong. In a comprehensive survey of smectic critical 
behavior in pure liquid crystals. Garland and Nounesis 
found that the critical exponents a, 7, and v± ob- 
tained from calorimetry and x-ray scattering studies dis- 
played complex but systematic trends as a function of 
Rm nil- Figures fnTal-fnTcl display the values of the 
critical exponents v\\, v±, and 7 of various smectic liquid 
crystals as a function of Rm, as originally assembled by 
Garland and Nounesis. Liquid crystals with small Rm 
generally have critical exponents that approach the 3D 
XY values; however, the values of the exponents do not 
change monotonically with decreasing Rm- The criti- 
cal exponents for pure 8CB {Rm = 0.977) are shown by 
the solid circles in the figures. Also in Figs. fTlTa)- fTlT c') 
are the effective critical exponents for 8CB confined to 
aligned aerosil gels, which have been placed on the fig- 
ures by assigning the smectics with quenched disorder an 
effective McMillan ratio. 



R"// = 0.977 - 



0.72p, 



(9) 



As the figures illustrate, this mapping between disor- 
der strength, as parameterized by aerosil density, and 
nematic-smectic coupling, as parameterized by McMil- 
lan ratio, places the effective critical exponents for 8CB 
confined to aligned aerosil densities roughly in line with 
trends of pure liquid crystals of varying Rm- We note a 
very similar linear mapping of aerosil density to effective 
McMillan ratio was introduced previously for 8CB and 
40.8 confined to isotropic aerosil gels [12, Ull- In these 
cases, the mapping was shown to collapse successfully 
effective heat capacity and order parameter exponents, 
a and /?, for various aerosil densities onto the trends of 
pure liquid crystals. The linear coefficient that optimized 
the mapping for the smectics confined to isotropic aerosil 
gels was smaller than that in Eq. 0.47 as opposed to 
0.72, suggesting that aligned aerosil gels are more effec- 
tive than isotropic gels in suppressing nematic-smectic or- 
der parameter coupling. Nevertheless, these linear map- 
pings provide a quantitative measure of how the presence 
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of the quenched disorder influences the smectic critical 
behavior. The results in Figs. [TITa^-fTlTc) thus extend 
the notion of an effective McMillan ratio to the effec- 
tive critical exponents v\f^, v^, and 7, making compre- 
hensive the analogy between disorder strength and R'^J/ 
for describing the behavior of the pre-transitional critical 
fluctuations of smectic confined to aerosil gels. Further 
theoretical efforts to understand the effects of quenched 
disorder could hence be a fruitful avenue for unraveling 
the outstanding mysteries surrounding the critical behav- 
ior of pure smectics. 
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